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Abstract 



^: 

' In D - dimensional dilaton gravitational model with the central charge deficit A 

Q> . the generalized Friedmann- type cosmological solutions (spatially homogeneous and 

isotropic) are obtained and classified. 

I— I ! 

r-- ■ Introduction 



One of the consequences of the string theory in its low-energy limit is the scalar dilaton 
field 0(x) [|l], in the space-time with the (critical) dimension D. The dilaton was involved 
Q\ I in various gravitational models. In particular, it makes non-trivial the 1 + 1 - dimensional 



D , authors [^] - searched 3 + 1 - dimensional solutions with the Friedmann- Robertson- 
^ ' Walker metric keeping in mind some form of compactification of other D — 4 dimensions. 

The dilaton field was recognized as a mechanism for the inflation EM that was required 



(2D) gravity that has fruitfully progressed during the last ten years P, ||, ^. 

On the other hand the dilaton field and other lowest-order corrections to Einstein gravity 
from the string theory were intensively used in cosmology - Some of the mentioned 



^ ■ for solving some cosmological problems. In Refs. |T3|, |T6|, |T9|, ^ the dilaton was used in 
. attempts to avoid the Big Bang singularity. 



A. Tseytlin |Tl[] obtained a set of cosmological solutions for the maximally symmetric 



space of an arbitrary dimension D with the dilaton, the central charge and the covariantly 
constant antisymmetric tensor field (for D = 3 or 4). 

The latter approach is used in the present paper. We consider and classify cosmological 
solutions in D- dimensional dilaton gravity with the central charge deficit A for the case of 
an arbitrary D. We use the D - dimensional generalization of FRW metric corresponding to 
homogeneous and isotropic spatial part with the constant curvature {k = 0, ±1). AH spatial 
dimensions are equal as it may be before the compactification. The purpose of our studying 
is to order the known 2D and higher D dilaton cosmological solutions (with adding new 
ones) and to classify them exhaustively. 

In this paper we don't consider the terms in the action with the antisymmetric tensor 
fields so the condition of the maximal spatial symmetry restricts the possible dimensions 



within D = 3 and D = A. In these cases the problem was studied by many authors [g, |TT 

0' S 
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1. Model and equations 

The dilaton gravitational model with the following action pi] is considered: 



S = J y^e-^-^^^) [R + 4(V0)^ + a] r+^x. (1) 



Here is the scalar dilaton field in ra+l- dimensional {D = n+1) pseudoriemann manifold 
Mi „ with coordinates x'^ and the metric tensor g^^, g = det \g^u\, R is the scalar curvature, 
(V0)2 = g^'-d^ct^d.ct). 

We use here and below the string frame for the action (|I]) keeping in mind that it may 



be conformally transformed g^j^y —>■ e '^'^g^Mui r = 4/(D — 2) to the Einstein frame [|I|, [IT], Q 
Equations of evolution for this model 

R^y -\{R + k) g^y + 20.^, = g^, [20;^ - 2(V0)2 

i(i?+A) + 0:^-(v0)2 = o 

are obtained from action (|l]) by the variation of the gravitational field g^y and the dilaton 
field accordingly. Here R^^ is the Ricci tensor, 0.^,^ is the covariant derivative. 
Transform the latter equations to the equivalent system 

R^y + 20.^,, = 0, , . 

i? + 4(V0)2 = A. ^' 

We search cosmological Friedmann's type solutions of Eqs. (^. It means that Mi^„ = 
RxS^ where S]^ is the homogeneous and isotropic manifold described by the scale factor a{t) 



Tl| . The manifold 5*^ has one of the following forms: n - dimensional sphere S"", pseudosphere 
S'!!:^, fiat space Sq = R^. In the spherical {k = 1) case a is the radius of S^. The metric in 
Mi^n = -R X "S"^ niay be taken in the form p5[ 



ds^ = de - a\t)[{dx^f + cl{x^){dx^f + ... + cl{x^) ■ cl{x^) . . . ■ cl{x''-^){dx'')% (3) 

{cosx, k = 1, 
1, = 0, 

coshx, k = —1. 

We substitute the metric (0) and = (j)(t) into Eqs. (Q) taking into account the following 
expressions (for all k): 

R^o = -nad/a"^, R\ = Rl = . . . = Rl = - [da + {n - l){d^ + k)] /a^, 
R = —n [2ad + {n — l)(a^ + k)]/a'^, d = da/dt. 

It reduces system (g) to three nontrivial equations 

-na/a + 20 = 0, (4) 
-a - (n - l)(a2 + A;)/a + 2d0 = 0, (5) 
n{n - l)(a^ + k)/a^ + 40^ - 4ra0d/a = A. (6) 

These three equations are not independent. Eq. is the "zero energy" constraint that 
is conserved as the consequence of Eqs. (^) and (||) and hence gives only a restriction on the 



initial values of a, d, 0, [O . 
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One may express from Eq. (^) 



a n — 1 2 ,\ 

- + —(a + k) 

a aa 



(7) 



and substitute it into Eq. (^. The resulting second-order equation (it is quadratic one with 
respect to a) may be transformed through the substitution 



a = q{a) 



to the following form: 



dq q . n — 1 q^ 
— k ± \ n— 



da a 



kn 



n 



+ A. 



aq 



(8) 



(9) 



To classify all solutions of the latter equation and system 
separately. 



we consider some cases 



2. Flat universe: k = 

In this case Eq. (P) is homogeneous one. It may be integrated for all values of A: 

_ a C^a^^^ - A/n 
It results in the following evolution of the scale factor: 

±l/v^ 

I tanni^^VArj 

a(t) = ao 



T 



tan(iA/w\r) 



A > 0, 
A = 0, 

A < 0. 



(10) 



Here C or its equivalent ao (oq = nC^/|A| for A 7^ 0) is the constant of integration in 
(|^). We use in Eq. (|TUp and below the time parameter 



T = ±it- to 



\t-t 



Oh 



because all solutions of system (P - @ are invariant with respect to time translation {to 
is the constant of integration) and time reflection. The signs ± in Eqs. (|11]) and ( |T0| ) (the 
latter is connected with the sign in Eq. (§)) are independent. 

The expression for (p corresponding to (^) results from Eq. (|^) 



20-200 



T 



sinh(|A/A 



r 



cosh(i-\/A ' 



sin(iv^w\r) 



cos(^-\/— Ar) 



-iTVn 



A > 0, 
A = 0, 

, A < 0. 



(12) 



The functions a(t) (|T^) and 0(t) (|T2D are exact solutions not only of Eqs. (|) and (||) but 
also of Eq. (^. This situation is more general. It is the non-trivial fact connected with the 
constraint character of Eq. (||) that all solutions of Eqs. (^ for all k, A and n satisfy Eq. (^). 
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It is necessary to supplement the family (0), (0) by the linear dilaton solution ^ 



a = const, 



±-VAt, A>0, k = 0, 



(13) 



case 



16, 20, Pl and were 



T^. They coincide with the solutions 



that was lost because of dividing on d in Eq. (|^). 

The solutions ([ToD , (|12D are well known in the AD or n 
considered for A > and an arbitrary D in Ref. 
obtained by Mueller for the manifold R x in the case when all radii of the n-dimensional 
torus are equal. The solutions of Veneziano ^ for the flat manifold with different scale 
factors generalize Eqs. (]T0|), (|T^). 

One can see that qualitative picture of solutions (|T0|) , (|12D and the structure of branches 
are similar for all n > 1. In the cases A > for all n we have two different branches of 
solutions (or 4 branches if we take two symmetric with respect to the change "+" to "— " in 
([TT| ) solutions as different ones [jT^)^. And there are only one branch for A < due to the 
symmetry properties of ([TT|) . 

The upper sign in Eqs. (|T0|) and (p!2| ) corresponds to an expansion (with increasing r) 
of the fiat dilatonic universe from the singular point r = with growing up from —oo, 
the opposite sign describes an evolution without the singularity a = but with a ^ oo, 
0— >-+ooatr— >-0. The value A = always corresponds to the power- law expressions for 
a and generalizing ones for the case n = 3. The presence of A > results in slowing 
down the above mentioned expansion and stabilization of the scale factor: a — oq = const, 
— > — oo at r — > oo. In the case A < we, v. v. have acceleration of evolution that is 
finished during finite time. 

Thus the described above role of A-term in the dilaton cosmology is in some sense opposite 
jTT| to that in usual Friedmann type n + 1- dimensional cosmology. In the latter case the 
evolution of the fiat universe [k = 0) with dust matter^ and with Einstein's gravitational 



Lagrangian J\g\ {R + 2 A) has the form 



a = const 



sinh I 
sin I 



riA 



2{n-l) 



2/n 



-nA 
2(n-l) 



T 



2/n 



A > 0, 
A = 0, 
A < 0, 



and A > accelerates the expansion. 

Note that for the 2D dilaton model (n = 1) all "cosmological" solutions are exhausted 
by Eqs. (|T0|) - (|1^) . It is naturally so the 1 - dimensional space has no an intrinsic curvature, 
hence k vanishes in Eqs. (^, (§) if n = 1. In this case the n = 1 solution (p!oD, (|T^) coincides 
with the D = 2 string "black hole" solution |T^, (taking the Euclidean form after the 
analytic continuation r 



ir, A 



A|Tg). 



3. Curved space k = ztl and A = 

For any value of k all solutions of system (^) - (^) may be obtained by integrating Eq. (|^). 
It was mentioned above that Eq. (^) in this case is satisfied too. But the most convenient 

"'^ Below for counting a number of the branches we don't differ two solutions connected only by the time 
reversal (11 ). Similarly, describing an expansion we always keep in mind that the corresponding solution with 
the contraction exists. If one differ these solutions he should multiply any mentioned number of branches by 
two. 

^In action (|l|) all matter is in "exotic" (dilaton and A-tcrm) form. 
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way to classify various types of the solutions for A = is to exclude the scale factor a from 
Eq. For this purpose we take the linear combination of these equations with the 

coefficients —1, n/a and 1. The resulting expression 



(with A = 0) is integrable: 



a 



(14) 



(15) 



Here C is an arbitrary constant. 

Before using the expressions (|14D or (|T5D we are to consider all solutions of system (H) - 
(|^) with = 0. For all k, A, n > 1 we have only two solutions of this type: the trivial case 
of ( |13|) a = Gq, (j) 



-- 00, = 0, A = and the solution with linear grows of a |]TT 
a = ao±t, = 00) k = —1, A = 0. 
Substituting Eq. ([T^) or ( 1T5|) into Eq. (^) we obtain the differential equation 

d 



n 



dt , 



(16) 



(17) 



Introducing the notation 



P = 0/0^ 

we integrate Eq. (|T7p in the following form: 



dt Q 



p + 



v^-1 



2(Vn-l) 



P 



2(v^+l) 



Here Q is the constant of integration that can take various values in three intervals resulting 
from the division of the axis — oo < p < oo by the singular points 



Pi 



P2 



A/n - r + 1 

Taking into account that pdt = —d (0)~^ we obtain the relation 



(19) 



dt 



'p — 1 2(^-1) \p — ^(^+^) sign (p — pi) sign [p — P2) dp. 



n — 1 



(20) 



It expresses t through the parameter p in quadratures in each of the intervals (— oo,pi), 

(Pl,P2), (p2,+00). 

Eqs. (|15D, ([18|) and (|3) let us express and a through p: 



n 



= 00 + ^ In 



p-pi 



P-P2 



(21) 
(22) 



Here the constants Ca, 0o and C in Eq. (|T3|) are connected: Ca = \CQe 



2(/.o|l/n 
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Expressions 



^ ) satisfy Eqs. (|^) - (|^) only under the following condition: 

{nQ/Caf = -ksign{p-pi) ■ sign{p - P2). (23) 



It means that the solution pOj ) - for the spherical case of A; = 1 exist only if p G (^1,^2) 
and for the case = — 1 it takes place if p G (— oo,pi) U {p2, +00). 

So for A; = 1, A = there is only one branch of solutions (we count the branches in 
accordance with the footnote in Sect. 2) represented in Fig. 1 as the evolution of the scale 
factor a in time (0) for various n = D — 1 (solid lines). The graphs for = 0(r) with 0o = 3 
are drawn by the dashed lines for the cases n = 2 and n = 9. For any n these solutions 
describe the expanding universe (a grows up from the singular value to +00 and rolls up 



from —00 to +00 when p changes between the points (O) pi and P2) during the finite time 



P2 



pi 



n 



V 



2(v/^-l)^ ^2(7^+1 



n 



'r 



n 



1 



(24) 



Here r(a;) is the Euler gamma function. In Fig. 1 the graphs a = a(r) and (p = 4'i''') with 
the equal values T = 1 for different n are drawn. Note that the solutions (EUp - (]2T|) for any 
fixed n > 1 and k = ±1 constitute the two-parameter family: the first parameter to in (p!T|) 
is trivial (a translation in t), the second is Ca, or Q, or T (they are connected by Eqs. (^) 
and (p^). The third parameter 0o contains Eq. (^2]). 
In the pseudospherical case k 



-1 the solutions (^OD - (p2|) are divided into two branches 
which are shown in Fig. 2. The first one (Fig. 2a) describes the expansion of a and growing 
0: < a < +00, —00 < (j) < (po, p < pi. The second branch with p > P2 (Fig. 2b) 
corresponds to evolution of a without the singularity a = 0: the value a decreases from 
+00 to the minimum a^m = Ca[\^/npl\P'-{^/np2y']~^^^ and then it grows up to +00 when 
(j) monotonically decreases from +00 to 0o- In these both cases the asymptotic t ^ 00 
behavior is — > 0o = const, a const + r, i.e. in this limit all solutions (pO]) - (^) for A; = 1 
tend to the linear evolution (|TBp. 

In the opposite (singular) limit p — *■ pi or p — > p2 for k = ±1 and for all branches we 
have an asymptotic dependence of the following two types: 



a ~ const ■ r 
I ^ const ■ r~^/^ 



const ■ T 
const ■ r 



r 
r 



+0, 
+0. 



(25) 
(26) 



Eq. ( P3| ) corresponds to the case p ^ pi and a that takes place for the branches 
represented in Figs. 1 (A; = 1) and 2a (A; = — 1). The case (p6D is realized for the solutions 
in Fig. 2b and in the limit t — T in Fig. 1. 

We see that the asymptotic behavior of the solutions (^) - (|2^ ) (A; = ±1, A = 0) and (0) - 
([T^ ) (A; = 0) is identical. So one can easily set the correspondence between the equivalent 
branches with the same limits (^) or (p^). 

There is only one^ dimension n = 4 (or D = 5) where an explicit solution (^OD - (^) 
exists: 



a{t) 



-T3 



r 



1/2 



0(t) 



T3 



(27) 



•^For the critical dimensions D — 10 or D = 26 the exponents in (| 
integrable through any Euler's substitution. 



are rational but this expression isn't 
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Here T = 8Q/\/3 in accordance with p^), r = ±{t — to) may be negative (unhke (]Tl|)). 

Expression ( ^7]) is rather simple and may illustrate the described above solution (|^) - 
(P^). The constant to in Eq. ( P?] ) is chosen so that the point r = is the image of p = p2 = 2/3 
and r = T is the image of p = pi = —2 for the continuous mapping the p-axis onto r-axis. 
So the points r = and t = T divide the r-axis into 3 parts, and the interval < r < T 
corresponds to the case k = 1. 

4. Curved space /c = ±1 and A 7^ 

In these cases we study various types of evolution a{t), (pit) by solving the system (H), 
(H) that may be rewritten as 



a = q, q = a ^ 



q — k{n — 1) ± qyn[q^ — k{n — 1)] + Xa? 



(28) 



To consider solutions of this system we use in this Section, in particular, numerical 
analysis so we have explicit solutions only in special cases. 

For example, for the dimensionality n = 4 or D = 5 (specified in the previous Section) 
and A > we obtained the exact particular solutions with the power-law evolution of the 
scale factor 

a = aoVr, e^-^-^*" = r ■ e'^^", = ^36/A, n = 4, k = ±1, A > 0. (29) 
For other values n the analogs of the solutions (^) also exist but the expressions 
a ^ aov^, e^'^-'^" ~ ^(«-2)/2 . ^k^Xr ^ ^2 ^ ^(n - 1)/VA, k = ±1, A > 0. (30) 



are only asymptotic at r ^ 00 [[TT]. In the opposite limit r ^ these solution have the 
asymptotic form (|25| ) for A; = ±1 or ( [26| ) for k = —1. 

We'll see below that the behavior of the mentioned solutions (|29|) and (|30D in the t — 00 
limit radically differs from that of all other solutions in the positive curvature case k = 1, 
A > 0. Otherwise, in the case /c = — 1, A > all solutions have the asymptotic form (|30D . 

Analyzing the solutions of the system (^) we are to note that there is the domain 

nq^ + Aa"^ - kn{n - 1) < 0, (31) 



in the phase ag-plane where the radicand in ( pq ) or (^) is negative. Hence, the graphs of the 
solutions in phase plane can not pass through the domain (|3ll) and cross its border — the 
2-nd order curve. 

For k = 1 and an arbitrary fixed value A the graph of an every solution of Eqs. (|28|) in the 
phase ag-plane contacts with the mentioned border of the domain (|31|), that is this border 
(the ellipse if A > or the hyperbola if A < 0) is the envelope of the graphs family. At the 
contact point (we denote its coordinates Oq, go) for each solution the sign before the radical 
in Eq. (pHj) changes from "— " to "+". So any graph of solution in the case k = 1 consists of 
two parts corresponding to these signs and smoothly connected at the point (ao,go). 

In Fig. 3 the graphs of the considered solutions with = 1 in the to-plane are represented 
for the case n = 3 and various A (Fig. 3a), and for the case n = 4 with fixed value A = 1 
(Fig. 3b). The parameter to in Eq- (O) is chosen so the value a = corresponds to t = 0. 
Under these conditions the curves in Fig. 3b are the one parameter family, and one may 
choose the value ao (if go > 0) as the parameter numerating the curves in the family. The 
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values Oo are shown near corresponding lines in Fig. 3b. The curves marked by symbols "4" 
and "5" in Fig. 3b correspond to the values go < at the point of sign reversal in Eq. (|28|) ; 
these numbers are maximal values of the scale factor a for these solutions. 

In Fig. 3a for the case n = 3 the solutions a = a{t) with the same value aQ = 1 for various 
values A (marked near the graphs) are shown. For the curve with A = 10 the value a = 1 is 
the maximum of a. We see that presence of A < accelerates the expansion in comparison 
with the case A = (Fig. 1) and A > — decelerates. But with increasing A (without 
exceeding the certain limit) the expansion remains irreversible — the values a and grows 
up to +00 during the finite time T(ao, A). The asymptotic behavior a{t) at t — >• T takes the 
form ( p6D or 

a ^ const ■ (T - t)-^/v^, e^"^ ~ const ■ (T - t)-^-^, t-^T. (32) 

If A reaches some critical value (depending on oq) the expansion becomes unlimited in time 
with the power-law behavior (|30D at t —>■ oo. In Fig. 3a {n = 3, = 1) this critical value 
slightly exceeds A = 2.732. For larger A the type of evolution changes — the expansion in 
some time interval is replaced by compression that is finished in the finite time T = T(ao, A) 
with vanishing a. The value (j){t) during the time T grows from —oo up to some maximum 
and then decreases to — oo. In the limit t ^ T the asymptotic relation ( pSj) takes place. 

Fig. 3b illustrates that the critical power-law solution (it takes the form (pOj) for n = 4) 
may be obtained for fixed value A > by choosing the value qq. Any small deviation from 
the critical value results in an evolution within the framework of any of the two mentioned 
types with finite lifetime T. 

For the case of negative spatial curvature k = —1 and A < each graph of solution in the 
phase ag-plane contacts with the border of the forbidden domain ( ^T]) with the sign reversal 
in Eq. (|28| ) (similarly to the case k = 1). The family of these graphs a = a{t) for fixed n = 4, 
A = — 1 and various values oq is shown in Fig. 4a. The coordinates of the contact point 
numerate the curves as in Fig. 3b. 

All the solutions with A < (Fig. 4a) have finite lifetimes T(ao, A) and the behavior (|3^) . 
But their behavior at t — depends on oq and A. If oq does not exceed some critical value 

= acr(A) then a, — -|-oo at t in accordance with (p6D (similarly to the case A = 
in Fig. 2b). In the opposite case qq > acr this limit is a, e^*^ — > (^) (compare with Fig. 2a). 
These two cases are divided by the solution with the critical value ao = cicr (Fig- 4a) that is 
close to the linear expansion (|T^) a ~ t, ~ 0o for t <^T. 



In the case = — 1, A > there is no real forbidden domain ( [3T| ) but the solutions with 
different signs in Eq. (|28D are connected by the time reversal t —t, q —q. So we have 
one family of solutions with the infinite lifetime and a ~ ^/T (pSf ) asymptotic behavior at 
t ^ oo [|ll[], that is shown in Fig. 4b for A = 1 and n = 0. The values g = d at the point 



a = Oo = 1 numerate the curves. For values q < gcr('20)A) (here gcr — 1) the solutions have 



the t — > limit (pq) and for q > qcr — the limit (pSf ). The critical solution with q = q, 



icr 



(Fig. 4b) has a ^ t, (f) c::^ (j)o dependence for small t. 
Conclusion 

The Friedmann type cosmological solutions in D - dimensional dilaton gravity with the 
action (0) are investigated in this paper. Different values of spatial curvature k = 0, ±1 
and various values A are considered. For all cases the solutions of the system (H) - (||) 
(cosmological solutions) are described and classified. They include as the expressions ( p!OD - 
dH), (HD (particular cases or generalizations of them were obtained in Refs. |0] - Hill), 
as new analytic solutions (|2g) - (^ , (gT]) , (0). 
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To classify various types of solutions we represent them in the following table: 

Table 1. Types of cosmological solutions. 





A < 


A = 


A > 


k = 1 


^ oo 


^ oo 


^ oo 
^ a ~ aoV^ 
0^0 


k = 


^ oo 


oo ^ a ~ aoT~-^/^ 

a = ao, = 00 
^ a ^ aoT^/^ 


oo ^ a ~ ao 
a = ao, = ^VAt 
^ a ~ ao 


k = -l 


oo ^ oo 
a ~ r ^ oo 
^ oo 


oo ^ a ~ r 
a = r, = 00 
^ a ^ r 


oo ^ a ~ ao-\/r 
a r ^ d ao^/T 
— s> a ~ aoV^ 



Here the symbol "0" denotes the beginning or the end^ of an evolution with the asymptotic 
behavior (pSD (a, e^"^ 0). The symbol "oo" corresponds to the asymptotics a,(f) —>■ +oo 
of the type (^) or (^). Note that for A < a lifetime of the dilatonic universe is always 
finite — the beginning and the end of the evolution are marked in Table 1 by the mentioned 
symbols (for = — 1 we have the critical solution with close to linear a ~ r, ~ 0o beginning 
of the expansion and the finite lifetime). 

Other symbols describing the end of an evolution in Table 1 correspond to an infinite 
lifetime and power-law evolution within one of the following types: ( pOD (a ~ y/r), (111 
(A: = 0, A>0), ®, {k = -l,a^r). 

-1, A > all solutions have the asymptotics aov^ 



ni, but for k 



For k = 

fixed A > and to only one critical solution (dividing different types of solutions) has such 
a behavior. 

Note that the considered cosmological solutions have identical qualitative features for 
various dimensionalities > 2 (n = 1 or D = 2 is reduced to the case k = 0). In this sense 
the classification in Table 1 is universal for all D. 
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